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g ! ABSTRACT 

■ Cosmological linear perturbation theory predicts that the peculiar velocity 

|-t j ■ V{ x ) an d the matter overdensity 5{x) at a same point x are statistically 

OO ■ independent quantities, as log as the initial density fluctuations are random 

Gaussian distributed. However nonlinear gravitational effects might change 
the situation. Using framework of second-order perturbation theory and the 
Edgeworth expansion method, we study local density dependence of bulk 



> 



velocity dispersion that is coarse-grained at a weakly nonlinear scale. For a 
. typical CDM model, the first nonlinear correction of this constrained bulk 

velocity dispersion amounts to ~ 0.35 (Gaussian smoothing) at a weakly 
nonlinear scale with a very weak dependence on cosmological parameters. We 
also compare our analytical prediction with published numerical results given at 



Or 



nonlinear regimes. 



^ ■ Subject headings: cosmology: theory — large-scale structure of 
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5-h ' 1. Introduction 
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The peculiar velocity field is one of the most fundamental quantities to analyze the 
large-scale structure in the universe {e.g. Peebles 1980). It is considered to reflect dynamical 
nature of density fluctuations of gravitational matter. The peculiar velocity field is usually 
observed using astrophysical objects {e.g. galaxies), as determination of distances is crucial 
for measuring peculiar velocities (Dekel 1994, Strauss & Willick 1995). There is a possibility 
that statistical aspects of the velocity field traced by these objects and that traced by dark 
matter particles might be different. This difference is generally called "velocity bias" and 
its elucidation becomes highly important in observational cosmology {e.g. Cen & Ostriker 
1992, Narayanan, Berlin & Weinberg 1998, Kaufmann et al. 1999). 



Velocity bias is often discussed numerically with making "galaxy particles" in some 
effective manners. But here we discussed a more basic phenomenon. It is known that 
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statistics of the peculiar velocity field depend largely on local density contrast. For example, 
both the single particle and pairwise velocity dispersions of dark matter particles are known 
to be increasing function of local density (Kepner, Summers, & Strauss 1997, Strauss, 
Cen & Ostriker 1998, Narayanan et al. 1998). Analysis of pairwise velocity statistics 
is interesting from theoretical point of views, and also very important in observational 
cosmology (Peebles 1976, Davis & Peebles 1983, Zurek et al. 1994, Fisher et al. 1994, 
Sheth 1996, Diaferio & Geller 1996, Suto & Jing 1997, Seto & Yokoyama 1998a, 1999, Jing 
& Borner 1998, Juszkiewicz, Fisher & Szapudi 1998, Seto 1999b, Juszkiewicz, Springel & 
Durrer 1999). But we do not discuss it here and concentrate on velocity field characterized 
by single point which is simpler to analyze theoretically. 

Linear perturbation theory predicts that the peculiar velocity V(x) and the density 
contrast 5(x) at a given point x is statistically independent, as long as the initial density 
fluctuations are random Gaussian distributed. Namely, the joint probability distribution 
function P(V, 5) can be written in a form as Pi(V)P 2 (£)- 

It is not surprising that the peculiar velocity of each particle is largely affected by 
nonlinear gravitational effects and shows local density dependence described above. But 
what can we expect for the smoothed (bulk) velocity that is field coarse grained at some 
spatial scale R? Due to nonlinear mode couplings, the relation Pi(V)P2(S) valid for linear 
theory must be modified and bulk velocity dispersion must also depend on local density 
contrast defined at the same smoothing scale R (see Bernardeau 1992, Chodorowski & 
Lokas 1997, Bernardeau et al. 1999 for the velocity divergence field). 

However, Kepner, Summers & Strauss (1997) showed from cold-dark-matter 
(CDM) and hot-dark-matter (HDM) N-body simulations that at nonlinear scales 
(0.77/i~ 1 Mpc < R < 4.88/i -1 Mpc), such a local density dependence was not observed (see 
Figs. 2 (a) and 3(a) of their paper). This is an interesting contrast to the behavior of velocity 
field traced by each particle, as described before (Kepner et al. 1997, Narayanan et al. 
1998). 

In this article, we investigate local density dependence of smoothed (bulk) velocity 
dispersion using framework of second-order perturbation theory. We calculate the first-order 
nonlinear correction of the constrained velocity dispersion. Our target is weakly nonlinear 
scale and somewhat larger than scale analyzed by Kepner et al. (1997). Since current 
survey depth of the cosmic velocity field is highly limited, our constrained statistics might 
not be useful in observational cosmology at present (e.g. Seto & Yokoyama 1998b, see also 
Seto 1999a). Our interest in this article is theoretically motivated one about nonlinear 
gravitational dynamics. 
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As the peculiar velocity field is more weighted to large-scale fluctuations (smaller wave 
number k) than the density field, perturbative treatment of smoothed velocity field would 
be reasonable at weakly nonlinear scale. Actually, Bahcall, Gramann & Cen (1994) showed 
that smoothed unconstrained velocity dispersion in N-body simulations are well predicted 
by linear theory even at smoothing scale R = 3/i~ 1 Mpc (see their Table 1). Second-order 
analysis by Makino, Sasaki & Suto (1992) also gives consistent results to their simulations. 



2. Formulation 

First we define the (unsmoothed) density contrast field S(x) in terms of the mean 
density of the universe p and the local density field p(x) as 

5(x) = P{X) ~ P . (1) 
P 

Many theoretical predictions of the large-scale structure are based on continuous fields, but 
observations as well as numerical experiments (such as, N-body simulations) are usually 
sampled by points where point-like galaxies (or mass elements) exist. In comparison 
of theoretical predictions with actual observations or numerical experiments, smoothing 
operation becomes sometimes crucially important to remove sparseness of particles' system. 
This operation is also important to reduce strong nonlinear effects which are difficult to 
handle theoretically. Thus it is favorable to make theoretical predictions of the large-scale 
structure including smoothing operation. We can express the smoothed density contrast 
field 5 R (x) and the smoothed velocity field V R (x) with (spatially isotropic) filter W{x, R) 
as 

5 R (x) = J dx' 3 5(x')W(\x - x'\,R), V R (x) = J dx' 3 V(x')W(\x - x'\,R). (2) 

As we discuss only the smoothed fields in this article, we hereafter omit the suffix R which 
indicates smoothing radius. 

The velocity dispersion Sy(5) for points x with a given overdensity 5(x) = 5 is formally 
written as 

2 _ {V{xf5 Drc [5{x) - 5]} 

^ v{d) ~ <W*(*)-*]> ' ( } 

where 5i)rc(') is Dirac's delta function and brackets (•) represent to take ensemble average. 

We assume that the initial (linear) density fluctuations are isotropic random Gaussian. 
At the linear-order we have V(x) oc VA^ 1 5(x) and (V(x)5(x)) = due to isotropy of 
matter fluctuations. This means that 5 and V at a same point are statistically independent 
quantities, as a multivariate probability distribution function (hereafter PDF) of Gaussian 
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variables is completely decided by their covariance matrix (e.g. Bardeen et al. 1986). 
Thus the constrained velocity dispersion £y(£) does not depend on the density contrast 
5 at linear order. However, nonlinear mode couplings would change the situation. Let 
us examine weakly non-Gaussian effects on Hy(5). We can express the first nonlinear 
correction of £y(<5), using framework of the Edgeworth expansion method (Cramer 1946, 
Matsubara 1994, 1995, Juszkiewicz et al. 1995, Bernardeau & Kofman 1995). This method 
is an excellent tool to explore weakly nonlinear effects of the large-scale structure induced 
by gravity. 

When a field F is defined by weakly non-Gaussian variables {A^(x)} with vanishing 
means, we can expand the expectation value (F) as (see appendix A) 

(F(.4 t , • • , A m) ) — (F) G + ^ £ (A^A X ) C {gJ^ M ;) o + Oi^F), (4) 

where (-) G is the expectation value under the assumption that variables {A^(x)} are 
multivariate Gaussian distributed, characterized by their covariance matrix [A^AJ). The 
quantity (A fl A u A\) c is the third-order connected moment of variables {^j, and we have 
(A^A u Ax) c = (A^AyAx) at third-order. The variance a 2 = 0(A 2 ) is the order parameter of 
perturbative expansion around the Gaussian distribution, and we can regard a 2 = (5 2 ) in 
this article. The denominator of £y(5) in equation (3) is nothing but the one point PDF of 
density contrast 5. From equation (4) we obtain the famous perturbative formula as follows 
(u = 5 /a) 

(SnrcMx) -6\) = ^(l + + 0(. 2 )) , (5) 

(e.g. Juszkiewicz et al. 1995, Bernardeau & Kofman 1995) and this is the most simplified 
version of the Edgeworth expansion. Here the function H n (v) = (— l) n e u2 l 2 (d/ di>) n e~ u2 1 2 is 
n-th order Hermite polynomial, and S is a parameter of order unity and called skewness 
(Peebles 1980, Fry 1984, Goroff et al. 1986, see also Seto 1999c), 

Due to the nonlinear correction term proportional to Sa, points with high-a overdensity 
are more abundant than the linear prediction by a Gaussian distribution. 

Next the numerator of Xy(5) is given by 

V(x) 2 5 Drc [5(x) ~ $}) = ^v^== f 1 + SaH Q (U) + CoH x (v) + 0(a 2 )\ , (7) 
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where <jy = (v 2 '") is the unconstrained velocity dispersion. The parameter C = 0(1) is 
defined by 

In the studies of the large-scale structure, the Edgeworth expansion or the third-order 
moments have been mainly discussed for scalar fields, such as density field S(x) or velocity 
divergence field V ■ V(x) (e.g. Chodorowski & Lokas 1997). Here we present analytical 
study for couplings of 5(x) and V(x), but numerical investigation of our method is also 
important as well as interesting. From equations (5) and (7) we obtain the constrained 
velocity dispersion £y(<5) up to the first-order nonlinear correction as 



{V(x) 2 8 Drc [5(x) - 5}} 
(5 Drc [5(x) ~ 8]) 

a- 



'v 



l + ^^ + CaH 1 (u) + 0(a^ 



= o* (l + SaH f v) +C5- SaH 6 3{l/) + Q(a^ 

= a 2 v (l + C5 + 0(a 2 )). (9) 

Note that our result does not depend on the skewness parameter S. Nonlinear effects 

appear through the quantity C. 

Next let us evaluate non-Gausssianity induced by gravity, using higher-order 
perturbation theory. We perturbatively expand the density and velocity fields as 

5(x) = S^x) + 5 2 (x) + 5 3 (x) + •••, (10) 
V(x) = Vi(aO + V 2 (aO + V- 3 (aO + ..., (11) 

where 8\(x) and V\(x) are the linear modes, 62 (x) and V2(x) are the second-order modes, 
and so on. We solve the following three basic equations (continuity, Euler and Poisson 
equations) order by order (Peebles 1980) 

^5(x) + -V[V(x){l + 5(x)}] = 0, (12) 
at a 

^-V(x) + -[V(x)-V]V(x) + ^-V(x) + -V(t ) (x) = 0, (13) 
at a a a 

V 2 (j)(x)-A-Ka 2 p(t)5(x) = 0, (14) 
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where a represents the scale factor. In these equations we have omitted explicit time 
dependence of fields for notational simplicities. We only discuss quantities at a specific 
epoch and there would be no confusion. 

Fourier space representation is convenient to analyze the nonlinear mode couplings. We 
denote the unsmoothed linear Fourier mode by 5u n (k). Then 5i(x) and V 1 (x) are written 
in terms of Su n (k) and W{kR), the Fourier transform of the filter function W(\x\,R), as 

/dk f dk %k 

— z exp(tkx)5 lm (k)W(kR), Vi(a:) = Hf J p exp(ikx)8 lin (k)W(kR), 

(15) 

where H(= dlna/dt) is the Hubble parameter and /(= din D/ din a, D: linear growth rate 
of density fluctuation) is a function of cosmological parameters Q and A, and well fitted by 

/ ~ ft 06 + A/30, (16) 

in the ranges 0.05 < < 1.5 and < A < 1.5 (Martel 1991). 

We define the linear matter power spectrum P(k) by 

(S H n(k)S lin (l)) = (2n) 3 5 3 Drc (k + l)P(k). (17) 

Then the dispersions a 2 and <7y are given by the following simple integrals of P(k) up to 
required order to evaluate the first nonlinear effects of Hy(6), 

a 2 = J ^P{k)W{kRf + 0{a% (18) 

4 = ^fJ^^PikWikRf + Oi**), (19) 

In this article we only use the Gaussian filter defined by W(kR) = exp[— (kR) 2 /2]. 

As shown in equation (8), the first-order nonlinear correction of the constrained 
dispersion Sy(5) is characterized by the factor C. We need the second-order modes 
52{x) and V~2{x) to calculate the first nonvanishing contributions of (V(x) 2 5(x)) . These 
second-order modes are given with linear mode 8u n {k) as (Fry 1984, Goroff 1986) 

r dkdl 

S 2 (x) = J exp[i(fc + l)x}5 lin (k)5 lin (l)H 2S (k, l)W(R\k + 1\), (20) 



(2k) 

(2tt) 6 \k + l\ 2 
where kernels H28 and H 2 v are defined as follows 



V 2 (x) = Hf J ^ ±S eX p[i(fc + l)x]5 Hn (k)5 Hn (l)n 2 v(k, l)W(R\k + l\), (21) 



, s I, x k-l (k l\ I, Jk-l\ 2 , x 
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n 2V (k,i) 



L + 



kl 

~2kJ 



'k r 
i + k 



kl 

kl 



(23) 



The factors K and L depend very weakly on cosmological parameters Q and A, and are 
fitted as (Matsubara 1995) 



A 



K{Q,X) ~ -IT 1 / 30 -— 1--Alog 10 ft 



7 







L(fi, A) 



3^-11/200 ^_ 



^Alog 10 fi 
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in the ranges 0.1 < Q < 1 and < A < 1. In the followings we neglect these weak 
dependence and simply put 

*-*-?■ 

Thus we can write down the third-order moment {V ■ VS) in the following form 
(V-VS) = (V 1 -V 1 S 2 ) + 2(V 1 -V 2 S 1 ) + 0(a 6 ) 



(24) 
(25) 

(26) 
(27) 



2H 2 f 



dkdl 

J2nf 



P(k)P(l) 



k -\n 2 s(k,i) + 2 k - {k 



I) 



k 2 l 2 



k 2 \k + l\ 2 



n 2V (k,i) 



xW(kR)W(lR)W(\k + l\R) + 0(a 6 ). (28) 

Due to the rotational symmetry around the origin, we can simplify the six dimensional 
integral dkdl to three dimensional integral dkdldu. Here, — 1 < u < 1 is the cosine between 
two vectors k and I and given by u — k ■ I /hi. Then we obtain the first nonvanishing order 
of (V ■ VS) as 

k 2 l 2 dkdl 



(V-VS) = 2H'fJ_ i duJ^£pip(k)P(l)exp[-e-f 



klu] 



x 



u 



5 



u 



' k I 



-TT 3 + X T + T- 



kl 



7 2, 
k + lu 



I k 



... J 3 u 
k(k 2 + l 2 + 2klu) \ 7 + 2 



'k r 

I + k 



-u 



(29) 



Note that the parameter C does not depend on the normalization of power spectrum (see 
eqs. [18] [19] and [29]). Furthermore, the factors Hf cancel out between <jy and (V ■ VS) and 
cosmological parameters are irrelevant for the factor C in our treatment (K = L = 3/7). 
Finally, we comment that even though the constrained dispersion Sy(5) changes by OyCS 
from the unconstrained value ay, the shape of the one-point PDF of velocity field with a 
given S keeps Gaussian distribution at the same order of nonlinearity. We can easily confirm 
this by calculating the ratio 



(Sl rc {V{x) - V)S Drc {S{x) - SJ) 
(S Drc (S(x) - $)) 



(2-3- 1 ttS^( ( 5))3/2 



cxp 



V 2 



2-3-^(5) 



+ 0(a 2 ) 



(30) 
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The factor 3 1 in the right-hand side arises from the dimensionality of the velocity vector. 
First-order correction is completely absorbed to the velocity dispersion Sy(5). 

3. Results 

In this section we numerically evaluate the parameter C for various power spectra. We 
first examine pure power-law spectra P(k) given by (n > — 1) 

P(k) = Ak n . (31) 

In this case C does not depend on the smoothing radius R, and we can simply put R = 1. 
Then the dispersions a 2 and a v are given as 

where r(n) is the Gamma function. As for the nonlinear coupling (V ■ VS) = 
(Vi • Vi^) + 2 (Vi ■ V2S1), we can write down the first contribution (Vi ■ Vi^) explicitly 
in terms of Hypergeometric functions as in the case of skewness parameter S (Matsubara 
1994, Lokas et al 1995). However, using mathematica (Wolfram 1996), we confirm that the 
second term 2 (Vi ■ V 2 Si) cannot be expressed in a closed form and numerical integration 
is required. These two terms diverge in the limit n — > — 1 where velocity dispersion a v also 
diverges, but the factor C approaches in this limit. 

In Fig.l we plot C as a function of spectral index n in the range — 1 < n < 2. The 
correction C is a positive and increasing function of n. This means that the velocity 
dispersion of high density regions are larger than that of low density regions. We have 
C = 0.314 at the scale-invariant spectrum n — 1. 

Next we examine C for a more realistic power spectrum P(k). We use CDM transfer 
function given in Bardeen, Bond, Kaiser & Szalay (1986) and assume that the primordial 
spectral index is equal to I. Then P(k) can be written as 

2 

[1 + 3.89g + (16.1g) 2 + (5.46g) 3 + (6.71g) 4 ]- 1/2 , (33) 

where q = fc/[(r7i)Mpc -1 ]. T is the shape parameter of the CDM transfer function and 
recent observational analyses of galaxy clusterings support T = 0.2 ~ 0.3 (e.g. Tadros et al. 
1999, Dodelson & Gaztahaga 1999). In Fig.2 we plot C as a function of smoothing radius 
R in units of [(r7i) -1 Mpc]. For this model the factor C depends weakly on the smoothing 



P(k) = Ak 



ln(l + 2.34g) 
2.34g 
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radius R and we have C ~ 0.30 at a weakly nonlinear regime R ~ 10/i _1 Mpc. In the limit 
R — > oo, C converges to 0.314 which is the same value of C for the power-law model with 
n = 1 presented in Fig.l. This is reasonable as we have 



P(k) 

lim — - — = const, (34) 

k->0 k 



for CDM models analyzed here. 



Our results obtained so far are the velocity dispersion for points constrained by the 
matter density contrast 5. One might have interest in the velocity dispersion constrained 
by the galaxy density contrast 5 g . Here, let us assume deterministic but nonlinear biasing 
relation for the smoothed galaxy distribution S g (x) and the matter distribution 5(x) as 

S g (x) = M(aO + b 2 (5(x) 2 - a 2 ) + 0(a 3 ), (35) 

where b\ and b 2 are some constants (e.g. Fry & Gaztanaga 1993). In this case we can easily 
show that the velocity dispersion T 1 y(5 g ) 2 for points x with S g (x) = 5 g is given by 

Z v (5 g ) = 4(1 + CSg/h + 0(a 2 )), (36) 

where the factors C and ay are same as those appeared in £ y (<5) (eq.[9]). Thus Sy(5 s ) 
does not depend on the nonlinear coefficient b 2 . This is also apparent when we write down 
S(x) using S g (x) and then insert this solution to equation (9). The factor proportional to b 2 
is higher effects than analyzed here. Note that in equation (36), the linear bias parameter 
bi appears by itself not in the usual form (3 = Q os /bi, and the overdensity 5 g dependence 
becomes smaller for larger b\. 

Kepner et al. (1997) numerically investigated the mean magnitude (|V(cc)|) of 
smoothed bulk velocity for points with given overdensity 5. Following the fact commented 
in the last paragraph of section 2, we can easily calculate this magnitude and obtain 

following result (see appendix B) 



Hv(6) = \I^ 2 V (5)(1 + O(o»)) = ^a v [l + ?± + 0(* 2 )) . (37) 

For a typical CDM model, our analytical result predicts that the magnitude /j,v(S) is 
expected to change ~ 0.155, according to the local density contrast 5. If we constrain points 
using overdensity of galaxies instead of that of the gravitating matter, the combination C5 
is replaced by C5 g /b\ in the above equation. 

Numerical results of Kepner et al. (1997) were given for CDM and HDM models 
with os = 0.67 normalization. Here as is the linear rms density fluctuation in a sphere of 
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8/i Mpc radius. They calculated the smoothed density and velocity fields with smoothing 
radius R at 0.77/i -1 Mpc < R < 4.88/i~ 1 Mpc. Thus their results are quantities at nonlinear 
regimes. It is true that simple application of our perturbative formula to their results 
would not be valid. However, surprisingly enough, the quantity /iy(5) shows almost no 5 
dependence in the range < 5 < 30. Q 

If (and thus /iy(<5)) shows no (^-dependence at nonlinear scale and our 

second-order analysis is valid for weakly nonlinear regime, we confront an interesting 
possibility. Namely, with parameterization of overdensity by normalized value v = 5/a, 
velocity dispersion does not depend on v at linear and nonlinear regime, but depends on it 
at (intermediate) weakly nonlinear regime. Furthermore, we should notice that the velocity 
dispersion of dark-matter particles (without coarse graining) depends largely on 5 (Kepner 
et al. 1997, Narayanan et al. 1998) .0 

To make clear understanding of these transitions, we need to numerically investigate 
the constrained dispersion Xy(£) in detailed manner with various smoothing length R, from 
linear to nonlinear scales. Performance of second-order perturbation theory for the velocity 
vector is also worth studying. 

4. Summary 

It is commonly accepted that the large-scale structure observed today is formed by 
gravitational instability from small primordial density fluctuations (Peebles 1980). In 
this picture, the peculiar velocity and the density contrast are fundamental quantities to 
characterize inhomogeneities in the universe. Linear analysis of cosmological perturbation 
theory predicts that, as long as initial fluctuations are random Gaussian distributed, the 
one-point PDF of the velocity field V(x) is statistically independent of the local density 
contrast 6(x). This is an important aspect of cosmological perturbation. 

However nonlinear gravitational evolution changes the situation. Due to nonlinear 
mode-couplings, the peculiar velocity field is no longer statistically independent of the local 
density field. Here we have investigated bulk velocity dispersion (£ y (5)) as a function of 
the local density contrast and calculated its first nonlinear correction using framework of 
second-order perturbation theory. Our target has been set at weakly nonlinear regimes 
where perturbative treatment must be reasonable. At present, survey depth of velocity field 



x It seems that the function fJ.v(S) in their figures shows extremely weak dependence of 5 around S ~ 0. 
2 Definitions of velocity dispersion in these two papers are not identical. 
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is highly limited and our constrained statistics might not be directly useful for observational 
cosmology. However, we believe that our theoretical study is important to understand one 
interesting aspect of the cosmic velocity field peculiar to its nonlinear evolution. 

We have shown that the first nonlinear correction of £y(<5) is proportional to the local 
density and strongly depends on the matter power spectra, but weakly on the cosmological 
parameters Q and A. For typical CDM model with primordially scale-invariant fluctuations, 
this first-order correction is about 0.35. If we use overdensity of galaxies 5 g , this correction 
term becomes 0.35 g /bi (the factor b\ is defined in eq.[35]). This dependence might be used 
to constrain the linear bias parameter b± itself (not in the usual form Q os /bi) in future 
peculiar velocity surveys. We have also shown that the constrained one-point PDF of 
velocity field keeps the Gaussian shape up to second-order of perturbation. First nonlinear 
effects are completely absorbed to the velocity dispersion £y(5). 

Numerical results by Kepner et al. (1997) have been compared with our analytical 
results. Their results show almost no 5- dependence, contrary to ours. However spatial scale 
of their analysis (strongly nonlinear regime) is largely different from ours (weakly nonlinear 
regime). In the forthcoming paper, detailed numerical analysis will be presented for various 
smoothing lengths R and power spectra (see also Seto 2000). Validity of second-order 
analysis as well as the Edgeworth expansion method for velocity vector would be also 
investigated numerically. 

The author would like to thank J. Yokoyama for useful discussion and the referee R. 
Juszkiewicz for helpful comments. This work was partially supported by the Japanese 
Grant in Aid for Science Research Fund of the Ministry of Education, Science, Sports and 
Culture No. 3161. 
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A. Weakly Non-Gaussian Averages 

In this appendix, we derive expression (4) for weakly non-Gaussian variables 
{A^} (/x = 1, • • • ,n) with (A^) = (e.g. Matsubara 1994). The partition function Z^J^) for 
a multivariate probability distribution function P(A IJ ) is defined as 

Z{J,) = f° d n P(A„) exp (i ^ J.A,) . (Al) 

J — CO v ' 

According to the cumulant expansion theorem (Ma 1985), the function lnZ(J M ) is a 
generating function of connected moments (A^ ■ ■ ■ A^ N ) c . Therefore, taking the inverse 
Fourier transform of equation (Al), the probability distribution function P(A tl ) is written 
as 

/ °° (-l)N qn \ 

^(AJ = exp £^-f- ]T (^■-A llN ) )Pg(AJ, (A2) 

\N=3 iV! Mi,-,MiV OA^A - ••OA^ J 

where the function P G (A^) is the multivariate Gaussian probability distribution function 
determined by a (n x n) correlation matrix M MI/ = (A^A,) as 

W = = exp (-^MM- 1 )^) . (A3) 
^/(27r)«detM V 2 ^ / 

If we have relations (A^ ■ ■ ■ A m ) c = 0(a 2N ~ 2 ) as predicted by higher-order perturbation 
theory, equation (A2) is perturbatively expanded as 

1 d 3 
P(A„) = P G {A,) - - £ (A,A u A x ) c P G (A„) + 0(5 2 ). (A4) 

fl V \ 

Evaluating the ensemble average of a field F(A IJ ) with this perturbative expression and 
taking partial integrals, we obtain expansion (4). 

One might consider that formula (3) with Dirac's delta function is somewhat 
indirect. We can obtain same results using formula for probability distribution function 
P(V\5) = P(V,S)/P(S) and evaluating P(V,5) and P(S) with expression (A4). 



B. Derivation of 

Here we derive expression (37). As in the case of the constrained velocity dispersion 
£y(<5) (eq. [3]), the quantity /iy(<5) is formally defined as 

s »p[M . (B1) 

(8 Drc [8(x) - S\) 



The r.h.s. of this equation is written as 

.(\V\Sh re (V(x)-V)6 Dre (6(x)-S)) 



(*) = / 

J — c 



dV- 



(5 Drc (5(x) - 6)) 



With the perturbative expansion given in equation (30) the r.h.s. of this equation 
evaluated as 



1 



(2 • 3-%£^))3/2 



roc 




/ dV 


exp ^ 


I — oo 





V 2 



2- 3-^(5), 



+ 0(a 2 ) 



I VI 



It is straightforward to obtain expression (37) given in the main text as follows 

~8~ 



Hv(S) 



3tt 



^(l + O^)) 




^(l + ^ + 0(a 2 ) 
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Fig. 1. — The second-order correction C for power- law matter spectra with Gaussian 
smoothing. 
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Fig. 2. — The second-order correction C for the CDM spectrum of Bardeen et al. (1986). 
We plot the factor C as a function of the smoothing radius R in units of (/ir) _1 Mpc. 



